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Abstract 

Employing the Foldy-Wouthuysen transformation it is demonstrated straightforwardly that the first 
and second Chern numbers are equal to the coefficients of the 2+1 and 4+1 dimensional Chern-Simons 
actions which are generated by the massive Dirac fermions coupled to the Abelian gauge fields. A 
topological insulator model in 2 + 1 dimensions is discussed and by means of a dimensional reduction 
approach the 1 + 1 dimensional descendant of the 2 + 1 dimensional Chern-Simons theory is presented. 
Field strength of the Berry gauge field corresponding to the 4+1 dimensional Dirac theory is explicitly 
derived through the Foldy-Wouthuysen transformation. Acquainted with it the second Chern numbers 
are calculated for specific choices of the integration domain. A method is proposed to obtain 3 + 1 
and 2 + 1 dimensional descendants of the effective field theory of the 4 + 1 dimensional time reversal 
invariant topological insulator theory. Inspired by the spin Hall effect in graphene, a hypothetical 
model of the time reversal invariant spin Hall insulator in 3 + 1 dimensions is proposed. 
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1 Introduction 



2 + 1 dimensional Dirac theory which was supposed to be physically unrealizable turned out to be 
indispensable to reveal the main features of graphene. At the Dirac points, in the low energy and long 
wavelength limit, the electrons of graphene effectively satisfy the 2 + 1 dimensional Dirac equation 
where the velocity of light c is substituted by the effective velocity v F [1, 2]. This had enormous 
impact on condensed matter physics. Based on the 2 + 1 dimensional Dirac Hamiltonian Haldane[3] 
discovered that it is possible to construct a model where an integer quantum Hall effect results without 
an external magnetic field. In spite of the fact that a magnetic field is not needed, it is a time 
reversal breaking (TRB) theory. Kane and Mele[4] incorporated the spin degrees of freedom into the 
Haldane's construction and formulated the time reversal invariant (TRI) spin Hall effect in graphene. 
This model[4] paved the way to the theoretical prediction of the topological insulator phase in real 
materials [5] which was observed for the first time in [6]. Topological insulators are usually defined to 
be ordinary insulators in the bulk which possess conducting states moving at the edge or boundary 
surface[7, 8, 9]. They can be classified by a new topological invariant called Z^flO]. 

Topological invariants emerge already in the quantum Hall effect. The Hall conductivity is given by 
the first Chern number [11, 12]. Moreover, in 2 + 1 dimensions a topological gauge theory is generated 
by integrating out the massive Dirac fermion fields coupled to Abelian gauge fields in the related path 
integral[13, 14, 15]. It is described by the 2 + 1 dimensional Chern-Simons action whose coefficient 
is the winding number of the noninteracting massive fermion propagator which is equal to the first 
Chern number resulting from the Berry gauge field [16, 17] of the quantum Hall states. Thus the 
Hall current can be acquired from a topological field theory which manifestly violates time reversal 
symmetry [18, 19]. One can also derive the TRI spin Hall current of the model introduced in [4] by 
calculating the related first Chern numbers[20]. 

The 4 + 1 dimensional Chern-Simons action which can be generated by the massive fermions 
coupled to Abelian gauge fields, is manifestly TRI and Qi-Hughes-Zhang[21] designated it as the 
effective topological field theory of the fundamental TRI topological insulator in 4 + 1 dimensions. 
They demonstrated that for the band insulators which can be deformed adiabatically to a flat band 
model, the coefficient of the effective action is equal to the second Chern number given by the related 
non-Abelian Berry vector fields. We will prove the equivalence of the coefficients of the induced 
Chern-Simons actions with the Chern numbers in a straightforward manner by employing the Foldy- 
Wouthuysen transformation. Qi-Hughes-Zhang constructed the 3 + 1 and 2 + 1 dimensional descendant 
theories by dimensional reduction from the 4+1 dimensional action of the massive Dirac fields coupled 
to external gauge fields. We mainly follow the approach of [21], though we only deal with the continuous 
Dirac theory and propose a slightly modified method to introduce the descendant theories which 
permits us to derive explicitly the related physical objects like polarizations. Moreover, we posit a 
hypothetical model of TRI spin Hall effect in 3 + 1 dimensions which may be useful to understand 
some aspects of physically realizable three dimensional topological insulators [22, 23, 24]. 

In the next section, we introduce the Berry gauge fields corresponding to Dirac fermions through 
the Foldy-Wouthuysen transformation. Section 3 is devoted to demonstrate in a straightforward 
fashion that the coefficients of the the Chern-Simons actions generated by integrating out the massive 
Dirac spinor fields which are coupled to the Abelian gauge fields in the path integral are equal to 
the first and second Chern numbers given, respectively, by the Berry gauge fields in 2 + 1 and 4 + 1 
dimensions. We consider the 2 + 1 dimensional topological field theory of the integer quantum Hall 
effect in Section 4 and recall how to construct the TRI spin quantum Hall effect in graphene which 
is a model of 2 + 1 dimensional topological insulator. Then, we present the dimensional reduction to 
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1 + 1 dimensions by obtaining the one dimensional charge polarization explicitly. In Section 5, we 
consider the 4+1 dimensional Chern-Simons field theory which was shown to describe the fundamental 
topological insulator. We derive the field strengths of the related Berry gauge fields to obtain the 
second Chern number and study dimensional reduction to 3+1 dimensions. By imitating the approach 
of [4] we theorized a hypothetical model in 4 + 1 dimensions which yields a TRI spin Hall current in 
3 + 1 dimensions by means of the dimensional reduction. Slightly modifying the approach of [21] we 
proposed a dimensional reduction procedure to 2 + 1 dimensions which provides explicit forms of the 
gauge field components which take part in the descendant action. In the last section we discussed the 
results obtained. 

2 The Foldy-Wouthuysen transformation and the Berry gauge field 



We consider relativistic electrons of charge e > with a characteristic velocity like the velocity of 
light c or the effective velocity v F as in graphene. To retain the formulation general we work in units 
h = c = v F = 1, as well as e = 1 and recuperate them when needed. Thus, the free, massive electrons 
are described by the Dirac Hamiltonian 

H = a-k + f3m. (2.1) 

In this section, vectors are d-dimensional, like the momentum k whose components are denoted kj; I = 
1, • • • , d. The Hamiltonian (2.1) can be diagonalized as 

UHU ] = E/3, (2.2) 

where E is the total energy 

E = ^/k 2 +m 2 , (2.3) 
and U is the unitary Foldy-Wouthuysen transformation 

v _ PH + E 

V / 2E(E + m) ' 

Through the transformation U a pure gauge field can be introduced as 

The Berry gauge field A follows by projecting (2.4) onto the positive energy eigenstates of the Dirac 
Hamiltonian (2.1). One can be convinced that eliminating the negative energy states is equivalent to 
the adiabatic approximation by revoking its similarity to suppression of the interband interactions in 
molecular problems [2 5]. Thus, we define the Berry gauge field as 

A=I+A U I+, (2.5) 

where 1+ is the projection operator onto the positive energy subspace. This definition of the Berry 
gauge field is valid irrespective of the dimensions of the Hamiltonian (2.1). In order to derive A 
explicitly let us adopt the following 2 N x 2 N ; N = [|], dimensional realizations of a and j3 
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Here p and the unit matrix 1 are 2^ 1 x 2 N 1 dimensional. In representation (2.6) the gauge field 
(2.4) becomes 

where ajj = —%[ai,aj]. Therefore, the Berry gauge field (2.5) results to be 

Although the field strength of (2.7) vanishes because of being a pure gauge field, the Berry curvature 

is non-vanishing in general. 

When we deal with 2n + 1 dimensional space-time coordinates where n = 1,2- •■ , the Berry 
curvature (2.9) can be employed to define the Chern number which is the integrated Chern character, 
as [26] 

Nn = 71^7 / d2Hk e A/2-/2„tr {T hh ■ ■ ■ Ti 2n _ lhn } . (2.10) 

For 2 + 1 dimensional systems the Berry gauge field is Abelian, so that T a b = dAb/dk a — dA a /dkb, 
where a, b = 1, 2, and the first Chern number is 

Nl = ^j ^^ab- (2.H) 
In 4 + 1 dimensions one introduces the second Chern number as 

N 2 = / d A ke ijkl tT{FijF k i}, (2.12) 

where i,j, k, I = 1, 2, 3,4. 

3 Topological field theories and the Chern numbers 

Field theory of electrons interacting with the external Abelian gauge field A a is given by the Dirac 
Lagrangian density 

C(i/>,$,A) =4>[ 1 a ( Pa + A a )-m]i;, (3.1) 

where a = 0, 1 • • • d. By integrating out the fermionic degrees of freedom in the related path integral 
one formally gets the action of the external fields as 

S[A] = -ilndet[i-f a (d a - iA a ) - m]. (3.2) 

For d = 2n one of the terms which it gives rise to is 

T[A n+1 ] = J [d qi ] ■ ■ ■ [dq n+1 ]A ai (Qi) ■ ■ ■ A a " +1 (9„+i)7r ai ... Qn+1 (<?i • • ■ ?„+i). 
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denotes the integral over the related phase space. At the order of first loop 

/d 2n+1 k 
( 27r )2n+l tr { Gf ( fc ) Aa i k ~ <H) G ( k - 9l) ■ ■ ■ A "n+l + Qn+i, k)}, 

where G{k) is the one particle Green function of the free Dirac field and A Q is the photon vertex. 
T[j4™ +1 ] g enera t es the In + 1 dimensional Chern-Simons term 

S e /f m = C n j d + xe 1 2n+1 A ai d a2 A a3 • • • d a , 2n A a2n+1 , (3-3) 

which can be taken as the effective topological action in the low energy limit. In the weak field 
approximation the coefficient C n can be written as [27] 



e 



aiPi—a„P„a. 



n+l 



° n ~ (n + l)(2n + l)! 9{1)/?1 ' ' ' d ( n )^ a ^"°^+^ qi • • • Qn+i)\ qi =o, 
where d^ n ) a = d/dln- Substituting the photon vertex with 

X a (k,k) = -id a G-\k), 

where G~ l (k) is the inverse of G(k), one can express C n as the winding number of the fermion 
propagator G [27]: 

( _j\n+l a 1 ...a 2n +i f j2n+l£, 

Cn = \J +1){2n + 1)l J 7 ^^tr{[G W a Ql G W - 1 ]...[G W a a2ri+1 G(/,)- 1 ]}. (3.4) 

In the rest of this section we will demonstrate that for the 2 + 1 and 4 + 1 dimensional Dirac 
theories the topological invariant winding numbers of the fermion propagator (3.4) are equal to the 
Chern numbers of the Berry gauge fields (2.10), up to constant factors. Obviously, these results can 
be obtained by other means; however we think that the following derivations are quite straightforward 
and clear. 

By virtue of the Foldy-Wouthuysen transformation U, one can invert (2.2) to write the Dirac 
Hamiltonian as H = EU^ (3U, which is suitable to express it in a projector form. Indeed, for k a = (w, k) 
we can write the inverse of the propagator as 

G~' 1 (k) = w + (E + ie)(P- - P+), (3.5) 

where the projection operators P + and P_, satisfying 

Pi = P+, Pl = P-, P+ + P- = 1, (3.6) 

are given explicitly as 

P + = U^I+U, P- = tfl-U. (3.7) 

The operators i+ and L_ project, respectively, onto the positive and the negative energy states which 
are obtained in representation (2.6) as 

1 \ _ / 

o or I o i 
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Now, (3.5) can easily be inverted to obtain the Green function as 

G{k) = + 



w - (E + ie) w + (E + ie) 

In the sequel we will not explicitly write the positive, infinitesimal parameter e unless necessary. 
Derivative of G~ l with respect to ko = w is 

Moreover, owing to the projector form (3.5) and the energy- momentum relation (2.3), it satisfies 

dG~\k) _k: dP + 

~^kY~ - e {p - ~ p+) ~ 2E m;- (3 - 9) 

The following relations between the projection operators P+ and P_ 

p OP- = _dP ±p = OP- p _. p JP^ = _pJP± = ^p 
+ dkj dki dki ' dki dki dki + ' 

can easily be derived by inspecting their basic properties (3.6). 
3.1 Relation between C\ and Ni 

In 2 + 1 dimensions integration of the massive Dirac fermions in the related path integral with the 
Lagrangian density (3.1) leads to the effective topological action 

S 2 e + / f [A] =dj d 3 xe^A^d u A p , (3.10) 

where /j,,v,p = 0, 1, 2. The coefficient C\ is given by 
1 f d 2 kdni 

Cl = ~T2 e " UP J ^{[GifcJ^GW-^GW^GW- 1 ]^^^^)- 1 ]}. (3.11) 

Making use of (3.8) we can express (3.11) as 

1 f d 2 kdm 

Ci = --e ab J -^-^triG^GWdaGik^GimGik)- 1 }. 

When we write the integrand explicitly by employing (3.9), obviously the quadratic terms in k{ vanish, 
so that we get 

1 fd 2 kdw. r , P + , P- v />1E ,2. 



c > = -V* J W tr{( (^ + <^n^ )(4£ 8 » F+oa ' p + 

-2k a \P- - P+}Gd b P + - 2k b d a P + G[P. - P+])}. 
One can observe that k a and k b terms combine to vanish 

26 ab tr{[ I -^3 - ^-J^] [fc(A P + + k b d a P + ]} = 0, 
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even before performing the w integration. The remaining terms, after revoking the infinitesimal 
parameter e, become 

/d 2 kdw Pa. P_ 

(2-k) 6 (w - E - ie) 2 (w + E + xe) (w - E - ie)(w + E + le) 1 

Integration over w yields 

Ci = ^e ab J d 2 ktr{P + d a P + d h P + }. 

Making use of definitions (3.7), one can easily observe that 

e ab tv{P + d a P + d b P + } = e ab tr{(I+Ud a U^(I + Ud b U^) + i + d a ud b uU + } 
= e ab ti{I + d a Ud b U^I + }. 

On the other hand, plugging the field strength (2.9) of the Abelian Berry gauge field (2.5) into (2.11) 
leads to 

Nt = J d 2 ke ab tv{I + d a Ud b UU + }. 

Therefore, we conclude that 

iVi 

Cl = (3.12) 

3.2 Relation between C 2 and N 2 

When we deal with 4 + 1 dimensions, the effective topological action (3.3) becomes 

S 4 e + f )[A] =C 2 J d 5 xe ABCDE A A d B A c d D A E , 

where A, B, ■ ■ ■ = 0, ■ ■ ■ , 4. C 2 is given by (3.4) as 

C 2 = -1- f ^^Le^^tiiGdAG^GdBG^GdcG^GdDG^GdEG- 1 }. (3.13) 
3 x 5! J (27r) & 

Plugging (3.8) and (3.9) into (3.13) leads to 

d 4 kdw r r kj 



C 2 = ^ jkl j d ^t.{GG[^iP--P + )-2E dl P + ]G 



[|(P_ - P + ) - 2Ed j P + ]G[^(P- - P + ) - 2Ed k P + \G[^(P. - P + ) - 2Ed l P + }}. 

Because of the symmetry properties the terms depending on k at the second or higher order vanish. 
Thus, we need to consider only the following terms 

C 2 = ^Tstijki J d 4 kdwtr{2E 4 GGd l P + Gd J P + Gd k P + Gd l P + 

-E 2 (kiGG[P- - P + ]GdjP + Gd k P + GdiP + + kjGGdiP + G[P- - P + ]Gd k P+GdiP + 
+k k GGd i P+Gd J P+G[P^ - P+}GdiP+ + kiGGdiP + GdjP + GdkP + G[P- - P+])}. 
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Inspecting the symmetry properties one can easily observe that the second and the fifth terms are 
summed to give a vanishing contribution. Similarly, one can show that gathered together contributions 
of the third and the fourth terms vanish by making use of the equalities 

d,P + G(P_ - P + ) = - d,P + , (P_ - P + )Gd t P + = d t P + (J^-^-). 

Hence, we get 

C 2 = y^tijki j d A kdwEh?[GGd i P + Gd j P + Gd k P + Gd l P + \. 
After restoring e this can be expressed as 

C * = g^tijhl J d 4 kdwEHv[( {w _ E _ i£ $ w+E+ie y2 + {w _ E _ i£ y [w+E+i£? )d i P + d j P + d k P + d l P + ] . 
By performing the w integration we obtain 

1 I' d A k 

C 2 = ~-e m J j^tr{d i P + P-d j P + d k P + P-d l P + }. (3.14) 

In terms of the explicit forms of P + and P- given by (3.7) we can express the integrand of (3.14) as 
e ijkl tr{d i P + P-d j P + d k P + P_d l P + } = e ijkl ti{(U^I + diU + d^I + U)P_ (U^I+djU + djtfl+U) 

{uU + d k u + d k uU + u)P-{u^i + diU + diuU+u)} 

= e ijkl tx{I + d i UP-d j tfl + d k UP-d l Un + }. (3.15) 

Now, the Berry gauge field (2.5) is non-Abelian and its curvature Fij can be written as 

Tij = U + diUdjUU + + iI + UdiU ] I + UdjUU + -i^j 
= U + diUP-djUU + - U + d j UP-d i U ] I + . 

Inserting it into the definition of the second Chern number (2.12) and inspecting (3.14) and (3.15) one 
concludes that 

Generalizing this method to higher dimensions is straightforward. 

4 2 + 1 dimensional theory and dimensional reduction to 1 + 1 di- 
mensions 

One can observe that by employing (3.12) in (3.10) the effective topological action of external gauge 
fields coupled to massive Dirac electrons living in 2 + 1 dimensions becomes 

S lff = ^ / d 3 xe^A^A p . (4.1) 

To calculate the related first Chern number (2.11), let us choose the representation a = (a x ,a y ), 
where a a are the Pauli spin matrices. This corresponds to set p a = (1, —i) in (2.8). Thus, the Abelian 
Berry gauge field can be written as 



Aa ~2E(E + m y (42) 
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It yields the Berry curvature 

(4 - 3) 

We plug (4.3) into (2.11) and perform the change of variable by (2.3) to express the related first Chern 
number as 

m f dE . , 

where the domain of integration D will be specified according to the model considered. If it is required 
to treat the E > and E < domains on the same footing, we can deal with 

m f m dE m f°° dE 1 fA . 

N ' = -iL&~L& =h <4 - 5) 

4.1 A model for 2 + 1-dimensional topological insulator 

Before presenting the graphene model of [4], let us briefly recall the interconnection between the 
quantum Hall effect and the Chern-Simons action in 2 + 1 dimensions. For electrons moving on a 
surface in the presence of the external in-plane electric field E = (E x , E y ,0) and the perpendicular 
magnetic field B = (0, 0, B z ) the Hall current is given by 

ja = <r H e ab E b . (4.6) 

Ignoring the spin of electrons the Hall conductivity is a topological invariant[ll, 12]: 

e H = jNl (4.7) 

Here Aq is the first Chern number resulting from the field strength J-~ B of the Berry gauge field obtained 
from the single particle Bloch wave functions which are solutions of the Schrodinger equation in the 
presence of the external magnetic field B z , integrated over the states up to the Fermi level Ep as 



J (2^ ^) 



A field theoretic description is possible in terms of the Chern-Simons action (4.1) with definition (4.8). 
In fact, the current obtained from the topological field theory (4.1), 

i — — f 8 U A p 
Jii — 2^ tz ^p u > 

gives for E a = d a Ao — doA a the Hall current (4.6). It also leads to the charge density jo = a H B, where 
the induced magnetic field is B = d x A y — d y A x . Note that B would also be generated by the Hall 
current (4.6) through the current conservation condition d a j a = -dtjo- We would like to emphasize 
the fact that the field theory (4.1) is not aware of the external magnetic field B z . External magnetic 
field is responsible of creating the energy spectrum whose consequences are encoded in the calculation 
of the first Chern number (4.8). 

By employing the Berry gauge field derived from the Dirac equation (4.2), we can still get the 
Hall conductivity as in (4.7) by an appropriate choice of the domain of integration D in (4.4). This 
construction does not necessitate an external magnetic field. For the first time in [3], Haldane described 
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how to obtain the quantum Hall effect without a magnetic field (vanishing in the average) through a 
Dirac like theory. To calculate the Hall conductivity following from the Dirac equation we let all the 
negative energy levels be occupied up to the Fermi level E F = m in (4.4), so that 

e 2 / m f m dE\ e 2 , . 

In [4], Kane and Mele incorporated the spin of electrons into the Haldane model[3] and proposed the 
following Hamiltonian for graphene 

Hq = a x r z k x + a y ky + ma z T z s z , (4.10) 

which leads to a TRI spin current. The mass term is generated by a spin-orbit coupling. The Pauli 
spin matrices (T xyz act on the states of sublattices. The matrix r 2 = diag(l,— 1) denotes the Dirac 
points K, K' which should be interchanged under the time reversal transformation. The other Pauli 
matrix s z = diag(l, —1) describes the third component of the spin of electrons which should also be 
inverted under time reversal transformation. Thus the time reversal operator is given by T = UK 
where we can take U = T y s y and K takes the complex conjugation as well as maps k — > —k. Therefore 
(4.10) is TRI: 

THgT- 1 = Hq. 

The Abelian Berry gauge field obtained from the Hamiltonian (4.10) can be written as [20] 

A a = 2E{ E + m f a b k b l T S Zl (4.11) 

where 1 T is the unit matrix in the t z space. The corresponding field strength is 

Jia = diag(4,jl,4,jl) = -^3 W (4.12) 

The indices tl and ± label, respectively, the third component of the spin and r z . The spin current 
defined as 

•s -t , -t 4 4 
3 =3+ + 3-~ 3+ ~3±, 

leads to the spin Hall current 

fa = VsHtabEb- 

The difference of the related first Chern numbers 

«E=m 



2vr 



r / d 2 k + + 



2 2' K 2 2 
gives the spin Hall conductivity a SH as 



E=-oo 

= (s + s) = 2, (4.13) 



= £ AJVl = h (4 - 14) 



10 



4.2 Dimensional reduction to 1 + 1 dimensions 

We would like to discuss dimensional reduction from 2 + 1 to 1 + 1 dimensions by slightly modifying 
the procedure described in [21]. The dimensionally reduced theory can be defined through the 1 + 1 
dimensional Lagrangian density 

£i+i (tp, tp, A) = $ \i' (p r + A r ) + -y 2 ( y - m] ip, 

where r = t,x and the external field C y (t,x) is the reminiscent of the gauge field A y . We define 
( y = k y + £, where k y is a parameter which permits us to deal with one particle Green function of the 
2 + 1 dimensional theory to derive the effective action of the external fields as in Section 3. In fact, 
integrating out the spinor fields %p in the related path integral yields the effective action 

S eff = G iD(ky) J dxdt((x,t)e rs d r A s . 

The coefficient G\D{k y ) is required to satisfy 

j G 1D {ky)dky = N U (4.15) 

where the first Chern number N\ is given by (4.4). Instead of the Cartesian coordinates we prefer to 
work with the polar coordinates k, 6, where k x = kcosO, k y = ksinO. Similar to (4.15) we would like 
to introduce G(9) satisfying 

/ G(6)d6 = N u (4.16) 
Jo 

and define the (1 + l)-dimensional effective action as 

S 1+1 = G(6) J dxdt((x,t)e rs d r A s . (4.17) 

Although it can be deduced directly from definition (4.16), we can also obtain G(6) by writing the 
components of the Abelian Berry gauge field (4.11) in polar coordinates: 

k 

Ae = -2E(E + m y Ak = °- 



The Berry curvature remains the same 

Tk6 = \ 



8(kA e ) dA k 



m 

2EP' 



dk dd 

and allows us to calculate explicitly G(6) as 

G{9) = i- J kdkF ke 
m f 

Now, one can define the one dimensional charge polarization [28, 29] P{9) by 

dP(9) 



dE _ Ni 
E?~2^' 



m =G(0). (4.18) 
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Adopting the first Chern number calculated in (4.5), N\ = 1, we solve (4.18) by 

P(9) = ±. (4.19) 

The physical observable is not directly the charge polarization given by P(0) but the adiabatic change 
in P(6) along a loop, which is equal to 

AP = P(2tt) - P(0) = 1. 

The (1 + l)-dimensional action (4.17) becomes 

S 1+1 = ^J dxdtA r e rs d s ((x,t), (4.20) 

for N± = 1. Action (4.20) leads to the current 

jr = ~^^^rsdsCip^i^)i 

known as the Goldstone-Wilczek formula[30] and gives the charge 

1 f d({x,t) , 1 _ 

In fact, it corresponds to solitons on polyacetylene with charge Q = 1/2 for £ changing from to n 
and Q = 1/3 for A( = 2ir/3 as it was obtained in [30]. 

5 4+1 dimensional topological insulator and dimensional reduction 
to 3 + 1 and 2 + 1 dimensions 

The topological field theory 

S$)[A] = J^L J d 5 xe ABCDE A A d B A c d D A E , (5.1) 

is designated as the effective action of the 4 + 1 dimensional TRI topological insulators in [21]. It 
follows from (3.3) by making use of relation (3.16). To derive the related second Chern number N 2 
we deal with the 4 + 1 dimensional realization of the Dirac Hamiltonian (2.1) which is provided by 

Observing that pi = (iai,ia2,ios, —1), the non-Abelian Berry gauge fields can be obtained from (2.8) 
as 

_ a 3 k 2 - 0-2^3 ~ <?ih . _ -q-3/ci + aik 3 - a 2 k A 
^~ 2E(E + m) ' 2 " 2E(E + m) ' 1 j 

_ q- 2 fcl ~ 0-1&2 ~ ^3^4 . _ aih + q- 2 fc 2 + cr 3 fc 3 

2£(£ + m) ' 2E(E + m) ' ib j 
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By definition the Berry gauge field corresponding to the 4 + 1 dimensional Dirac Hamiltonian can also 
be derived by considering the explicit solutions of the Dirac equation as it was done in [31]. They 
work in the chiral representation, so that the Berry gauge field components which they obtain differ 
from (5.3), (5.4). 

One can show that the field strength components Tij = dAj/dhi — dAi/dkj — i[Ai,Aj], are 

[a 3 (-E(E + m) + k\ + k\) + a 2 (kik A - k 2 k 3 ) - a 1 (k 2 k 4 + hk 3 )] , 



?12 

•^23 
-^24 
•7"34 



2E 3 (E + 


m) 


1 




2E 3 (E + 


m) 


1 




2E 3 (E + 


m) 


1 




2E 3 (E + 


m) 


1 




2E 3 (E + 


m) 


1 




2E 3 {E + 


m) 



a 2 (E(E + m) — k\ — k 3 ) + a\(k\k 2 - k 3 k 4 ) + u^ikiki + k 2 k 3 )] , 
ai(E(E + m)-kl- k\) - o 2 (k x k 2 + k 3 k 4 ) - a 3 (k x k 3 - k 2 k 4 )] , 
o\{—E{E + m) + k\ + k\) - <J 2 (k\k 2 + k 3 k 4 ) - <J 3 (k\k 3 - k 2 k 4 )] , 
o 2 {E(E + m) - k 2 - k\) - a\(k\k 2 - k 3 k 4 ) - a 3 (k\k 4 + k 2 k 3 )] , 
a 3 (E(E + m) -k 2 - k\) + a 2 (k x k 4 - k 2 k 3 ) - o\(k 2 k 4 + k ± k 3 )] . 



Plugging them into (2.12) and taking the trace yield 



N 2 



47T 2 



m 

2~W 



)d 4 k. 



(5.5) 



To calculate it explicitly, we would like to deal with the four dimensional polar coordinates given by 
k\ = k cos <pi, k 2 = k sin <fii cos (j) 2 , k 3 = k sin <p\ sin <p 2 cos <fi 3 and k 4 = k sin <pi sin (p 2 sin <p 3 , where the 
angles 01,02,03, respectively, take values in the intervals [0,7r], [0,7r], [0,27r]. The volume element is 
d A k = k 3 sin 2 (pi sin 4> 2 dkd<pid(j) 2 d(j) 3 - Hence, after the change of variable by (2.3), one can show that 
(5.5) can be written as 

3m f m 2 — E 2 , ^ . 
N * = -} D —E^ dE - < 5 ' 6 > 
When D is taken to be an overlap of the E > and E < domains, we may deal with 



No 



3m f r 



m 2 -E 2 ^ 3m f°° m 2 — E 2 



(5.7) 



5.1 Dimensional reduction to 3 + 1 dimensions 



Dimensional reduction of the 4 + 1 dimensional effective action given by (3.1) to 3 + 1 dimensions can 
be described by the Lagrangian density 



C A+1 bP, ■>!>, A] = ij) 7 Q {p a + A a ) +j 4 e-m 



(5.8) 



where a = 0, • • • ,3. The external field 9{x a ) is the reminiscent of the gauge field A 4 . ip, ip fields can 
be integrated out as in Section 3 through the one particle Green function of 4 + 1 dimensional theory 
introducing the parameter k 4 by setting 9 = k 4 + 6(x a ). By keeping track of the phase space volume 
one can obtain the 3 + 1 dimensional effective action as 



b eff ~ 



G 3D (k 4 ) 
4tt 



J d 4 x6e al3 ^d a 



ApdjA^, 



(5.9) 
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where the coefficient is given through the condition 

G 3D (k 4 )dh = N 2 . (5.10) 



We would like to modify this construction by working with the four dimensional polar coordinates 
and proposing that the action describing the descendant theory is given by 



Ss+i = G^s) J ^ xee aMv daApdiA ^ 



whose coefficient, like (5.10), is required to satisfy the condition 

1^ G 3 (fa)dfa = N 2 . 
Jo 

Thus, the coefficient Gs(fa) can be obtained as 

G 3 (03) = ^2 J eijkitr^ijT^k 3 sin 2 fa sin fadkdfadfa = 7^, (5.11) 

with definition (5.6) of the second Chern number N 2 . 

Similar to the one-dimensional charge polarization (4.18) one can associate the coefficient G 3 {fa) 
to ^3(^3) through the relation [21] 



L 



;: :i r G3 (fa) d fa=N 2 . 



lo d<t>3 Jo 
Hence the "magnetoelectric polarization" can be obtained as 

Wa) = ^fa- (5.12) 

Observe that like the one-dimensional case, for A fa = 2ir it changes by AP3 = 1 if we choose N 2 = 1 
as it is calculated in (5.7). ^3(^3) depends linearly on fa due to the fact that the second Chern 
character corresponding to free Dirac particle depends only on k. Interacting Dirac particles may give 
rise to polarizations which would not be linearly dependent on fa. 

By inserting definition (5.11) into (5.9), the effective action becomes 

S 3+1 = ^J d^xOe^daApd^. (5.13) 

It can be written equivalently as 

= ^ J d 4 xP 3 (9)e a ^d a A^A v , 

where P 3 (0) = N 2 6/2ir. This describes the axion electrodynamics which is invariant under the shift 
9 -> 9 + 2-ir [32, 33]. 
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5.2 A hypothetical model for 3 + 1 dimensional topological insulators 

In spite of the fact that the underlying topological gauge theory (5.1) is manifestly TRI, the theory 
given by the descendant action (5.13) is TRB except for the values 6 = 0,7r. Nevertheless, we may 
deal with the TRB action (5.13) but introduce a TRI hypothetical model generalizing the spin Hall 
effect for graphene[4]. The current following from action (5.13) is 

f = A^a^. 

Assuming 9 = 9{z) and considering the in-plane electric field E a (x,y); a = 1,2, we obtain the 
current [34] 

N 2 
(2vr) 5 

The Hall current can be introduced by integrating (5.14) along the coordinate z as 



3a = 7K^2 d zO{z)e ah E h {x,y). (5.14) 



J a (x,y) = J j a dz = a H e ab E b (x,y). (5.15) 



It leads to the surface Hall conductivity a H [9, 21] 

a « = TW?! dz6dz = ^^Y M - (5 ' 16) 

Obviously we defined AO = 9(oo) — 9(—oo), which is non- vanishing for an adequate domain wall or at 
an interface plane between two samples. 

Now we should define the second Chern number (5.6) appropriately. We suppose that all negative 
energy states are occupied till the Fermi level taken as the first positive energy value E F = m, so that 
we get 

Considering a plane of interface which yields A9 = 2ir the Hall conductivity becomes 

e 2 

aH= 2h- 

In representation (5.2), the 4 + 1 dimensional Dirac Hamiltonian (2.1) is TRI where the time 
reversal operator can be taken as T 4+1 = a2Ct^K. However, the Hamiltonian corresponding to action 
(5.8) for A a = 0, 

H 3+1 = at ■ k + a 4 (9 + mfi, (5.18) 

violates time reversal symmetry. We will present a hypothetical model which is time reversal invariant 
emulating the spin Hall effect for graphene. Let on act on sublattices with two Dirac points. Assume 
that around these points which are interchanged under time reversal transformation, electrons are 
described by the Hamiltonians as in (5.18). Moreover, the third component of the spin given by the 
Pauli matrix s z = diag(l, —1) is included and conserved. Thus, we propose to consider the Hamiltonian 

H 3+1 = a ■ k + a.J) + T z s z /3m, (5.19) 
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where t z = diag(l, —1) and in terms of q« and j3 given by (5.2) we denned 

5j = (ai,T z a 2 ,a 3 ,a 4 ). 

Now, as in Section 4.1, the time reversal operator interchanging the Dirac points and the third compo- 
nents of the spin can be defined by T = T y s y K, so that (5.19) is TRI. Obviously, we can obtain (5.19) 
through the dimensional reduction from the 4+1 dimensional action corresponding to the following 
free Hamiltonian 

H = 5i- ki + T z s z Pm = Amg{H^ + ,H^-,H^ + ,H^). (5.20) 

As we show in Appendix, the four dimensional Hamiltonians defined by (5.20) correspond to the 
second Chern numbers 

jv| + = jv£ _ = -N^ + = -N^ = N 2 , 

where N 2 is given by (5.6). Repeating the procedure yielding (5.15)-(5.17) in the presence of a domain 
wall we can obtain the dissipationless spin current as 

J S a = 4 + + 4" " 4 + - 4~ = VsHe ab E b (x, y), 
with the spin Hall conductivity 

e / vt+ , Art- _ AT l+ _ atI-^ e 



for A9 = 2tt. It is equal to the spin Hall conductivity for graphene (4.14). 

5.3 Dimensional reduction to 2 + 1-dimensions 

The 2 + 1 dimensional Lagrangian density 

C 2+1 [ip,i>,A] = ip [f 4 (p M + A^) + 73 C 3 + 74C4 - m] ij>, 

describes the dimensionally reduced theory. The fields Q 3 {x p ), Q^x^) are the reminiscent of the gauge 
fields A3, A 4 , of the 4+1 dimensionally theory whose action is given by (3.1) for d = 4. By setting 
Cd,{xp) = k 3 + 4>{Xfj) and (, 4 {x p ) = k 4 + #(x^), where k 3 ,k 4 are parameters playing the role of the 
momentum components in one particle Green functions, one can follow the approach of Section 3 to 
derive the 2 + 1 dimensional effective action as 



S 2 e + f } = G 2D (k 3 ,k 4 ) J d z xe^ p A^d p l 



/' 



Its coefficient should fulfill the condition 

G 2D (k 3 ,k 4 )dk 3 dk 4 = N 2 . (5.21) 
As in Section 5.1, we consider the four dimensional polar coordinates and propose that the action 

S 2+1 = G 2 (4> 2 ,4>3) J d^e^A^dpO, (5.22) 
describes the 2 + 1 dimensional descendant theory. Obviously, like (5.21) we pose the condition 

fir c2tt 



JO JO 



d(f> 2 d<p 3 G 2 ((f) 2 ,(t)3) = N 2 . 
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This can be solved by 

N 2 

Proceeding as in [21] we introduce the vector field 

= ngd^e + n^d^, 

however by adopting the definitions 

N 2 ^2 

0,g = — — cos 4>, = — —9 sin </>. 
The field strength of the field fi^ is 

N 2 

d^9, v - d v 9.^ = — sin <j) {d u 9d^(j) - d^Od^) . 

Let us define <f> = 02 + an d = 03 + 9 as slowly varying fields, so that at the first order in derivatives 
we can write G 2 ((p,9)d ll 9d u (f) « G 2 ((f)2, <t>z)dn9d v <i>. Therefore, action (5.22) can be written in the form 

S 2+1 = ^J cPxe^A^Qp. (5.23) 

The current generated by the field Q a , a = 1, 2, is 

i° - 6 e J? 
J a - ^tab&b, 

where the electric field is given by E a = d a A$ — doA a . It can be interpreted as the spin current 
yielding the spin Hall conductivity a SH = e/2n. Hence, by attributing the adequate time reversal 
transformation properties to the gauge field Q^, action (5.23) corresponds to the TRI 2 + 1 dimensional 
model of [4] which we discussed in Section 4.1, 
Action (5.23) generates the electric current 

f = ^ UP£I ^p- 
For fields satisfying = <f>(x), 9 = 9{y) it gives the total charge 

jv 2 r n 2 r f 2w 

Q = e — / dxdy sin 4>d x 4>d v 9 = e — / sin0<i0 / d9 = eN 2 . 
4vr J 4n J J 

On the other hand the three dimensional Skyrmion field n coupled to Dirac fermion in 2 + 1 dimensions 
yields the current [35] 

j* = ^: W n • dUfl x dPn - 

The Skyrmion field configuration discussed in [36] satisfying n 2 = 1 possesses the charge Q T = 2e. 
Hence, if we deal with N 2 = 1, the Skyrmion theory can be described for the field configurations 
satisfying 

sin 4> (d u 9d^4> - d^9d v <j)) = -n ■ d v n x <9 M n, 
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which leads to = 2j fM . In principle, this condition can be solved to obtain n in terms of the fields 4> 
and 6. 

Observe also that for the field configurations <p = <t>(t), = 9{y) the net charge flow in x direction 

is 

J dtdyj x = -N 2 . (5.24) 

Moreover, we can introduce the magnetoelectric polarization in the form given (5.12) by defining it as 

r No 

Then the pumped charge (5.24) can also be written as AQ = f dP% = N 2 which gives AQ = 1/2 for 
N 2 = 1/2, as it is given for the Hall effect (5.17). 

6 Discussions 

The Foldy-Wouthuysen transformation which diagonalizes the Dirac Hamiltonian is proven to be a 
powerful tool to perform calculations in the effective field theory of the 4 + 1 dimensional TRI topo- 
logical insulator. The Foldy-Wouthuysen transformation is employed to obtain the Berry gauge fields 
of Dirac Hamiltonians and through them we derived the first and second Chern characters explicitly. 
Then we demonstrated in a transparent manner that the winding numbers of Dirac propagators are 
equal to the coefficients of the effective Chern-Simons actions in 2 + 1 as well as in 4 + 1 dimensions. 
This construction can be generalized to higher odd dimensions straightforwardly. In the line of the 
graphene model [4] we introduced a hypothetical model leading to a dissipationless spin current in 
3 + 1-dimensions. It can be helpful to understand some aspects of three dimensional TRI topological 
insulators if we can show that it is somehow related to some realistic models. Moreover, it seems that 
in terms of our explicit constructions one can discuss Z 2 topological classification of TRI insulators in 
a tractable fashion. In principle our approach can be generalized to interacting Dirac particles where 
the related Foldy-Wouthuysen transformation at least perturbatively exists. 



Appendix 

The Hamiltonian (5.20) which comprises t z and spin degrees of freedom denoted ± and ti, respectively, 
yields the 4+1 dimensional Dirac Hamiltonians 

= a±ki + a 2 k 2 + 03/03 + 04/04 + m/3, = a\k\ — a 2 k 2 + 03/03 + 04/04 — m/3, (A.l) 

= a±k\ + a 2 k 2 + 03/03 + 04/04 — m/3, = a\k\ — a 2 k 2 + 03/03 + 04/04 + m/3. (A. 2) 

Let us first consider the two spin up Hamiltonians (A.l). They yield slightly different non-Abelian 
Berry gauge fields 

4* = 2E{ E + m) (±^ - ^3 + ^ 4 ), 4* = 2E{ E + m ) iT(J3kl ± aih - a2h) > 

A ^ = 2E(E + m) {a2kl T aih T a3h) > A ^ = 2E(E + m) i±aikl + ^ ± a * h) - 
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The corresponding field strengths can be calculated as 

[^a 3 (E(E + m)-kf- kl) + a 2 (k 1 k 4 - k 2 k 3 ) T cri{k 2 k A + k x k 3 )\ , 
[a 2 (E(E + m)-kl- fcjj) ± <7i(fciA& - k 3 k 4 ) ± ^(fcifo + fofo)] > 
[±<7i(.E(.E + m) - kf - kl) - u 2 {kik 2 + k 3 k 4 ) T cr 3 (fcifc 3 - k 2 k 4 )] , 
[t<ti(-E(^ + m) - kj - kl) - a 2 (hk 2 + k 3 k A ) T <J 3 {kik 3 - k 2 k 4 )] , 
[cr2(£'(-E + m) — kl - k\) =F o"i(&i&2 - £3^4) T ^3(^1^4 + ^3)] , 
[±cr 3 (.E(.E + m) - fef - fc|) + o 2 (k x k± - k 2 k 3 ) T v\{k 2 k A + hk 3 )] . 
Although they are different, they generate the same second Chern number equal to (5.5): 



T-t± 

• r 12 ~~ 


1 


2E 3 (E + m 


7rt± 

^13 ~~ 


1 


2E 3 (E + m 


T-t± 

TIT = 


1 


2E 3 (E + m 


x-t± _ 
^23 ~~ 


1 


2£ 3 (£ + m 


•^24 — 


1 


2£ 3 (£ + m 


•^34 — 


1 


2E 3 (E + m 



N. 



32tt 2 / 



d A ke ijkl tx 



Trt± T-t± 



3 
4^2 



)d 4 £:. 



(A.3) 



The non-Abelian Berry gauge fields corresponding to the two spin down Hamiltonians (A. 2) can 
be shown to satisfy 

AfihMMM) = (-l^A^ik^k^k^-h), 

without summation over the repeated indices. Thus, the components of the related Berry curvature 
are 

rf(ki,k 2 ,k 3 ,k 4 ) = {-l) Sii+ ^^{k u k 2 ,k 3 ,-kA). 



They yield the same second Chern number which is given by (A.3) up to a minus sign: TV, 



2 



-N. 
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